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Abstract 

We consider the wetting transition on a planar surface in contact with a 
semi-infinite fluid. In the classical approach, the surface is assumed to be solid, 
and when interaction between solid and fluid is sufficiently short-range, the 
contribution of the fluid can be represented by a surface free energy with a 
density of the form $(^5), where ps is the limiting density of the fluid at the 
surface. 

In the present paper we propose a more precise representation of the surface 
energy that takes into account not only the value of ps but also the contribution 
from the whole density profile p{z) of the fiuid, where z is coordinate normal to 
the surface. 

The specific value of the functional of ps at the surface is expressed in mean- 
field approximation through the potentials of intermolecular interaction and 
some other parameters of the fluid and the solid wall. 

An extension to the case of fluid mixtures in contact with a solid surface is 
proposed. 



1. Introduction 

The phenomenon of surface wetting is a subject of many experiments. They 
have already been used to determine many important properties of the wetting 
behavior for liquids on low-energy solid surfaces ^. To gain a theoretical ex- 
planation of the phenomenon of wetting, a generalized van der Waals model is 
often used ^. 

In the 1950's Zisman ^ developed an experimental method of characterizing the 
free energy of a solid surface by the measurement of the contact angle with 
respect to the liquid-vapor surface tension and by changing the test liquid. 
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More recent measurements done by Li and all ^''"^ improved the understanding 
of this problem. While the contact angle and surface tension are macroscopic 
quantities, they have their origin in molecular interactions. 

In the following, we use a mean-field theory to investigate how the surface 
energy is related to molecular interactions. The approximation of mean field 
theory is too simple to be quantitatively accurate. However it does provide a 
qualitative understanding and allows one to calculate explicitely the magnitude 

of the coefficients in our model. 

In 1977, John Cahn^ gave simple illuminating arguments to describe the inter- 
action between solids and liquids. His model is based on a generalized van der 

Waals theory of fluids treated as attracting hard spheres^ . It entailed assigning 
to the solid surface an energy that was a functional of the liquid density " at the 
surface"; the particular form = —71 ps + — Ps of this energy, where 

ps is the fluid density at the solid wall, is now widely known in the literature 

and is due to Nakanishi and Fisher ^. It was thoroughly examined in a review 
paper by de Gennes ^ . To account for the wetting behavior of liquids on solid 
surfaces, one needs to know $(ps); a major object of the present paper is to 
obtain explicit formulas for the coefficients 71 and 72 , expressing them in terms 
of the parameters in assumed microscopic interaction potentials. 

Three hypotheses are implicit in Cahn's picture. 

i) For the liquid density to be taken to be a smooth function p{z) of the 
distance z from the solid surface, that surface is assumed to be fiat on the 
scale of molecular sizes and the correlation length is assumed to be greater than 
intermolecular distances (as is the case, for example, when the temperature T 
is not far from the critical temperature Tc). 

a) The forces between solid and liquid are of short range with respect to 
intermolecular distances. 

Hi) The fluid is considered in the framework of a mean-field theory. This 
means, in particular, that the free energy of the fluid is a classical so-called 

" gradient square functional" . 

The point of view that the fluid in the interfacial region may be treated as bulk 

phase with a local free-energy density and an additional contribution arising 
from the nonuniformity which may be approximated by a gradient expansion 
truncated at the second order is most likely to be successful and perhaps even 
quantitatively accurate near the critical point ^ . Some numerical calculations 
based on Cahn's model and comparison with experiments can be found in the 
literature 

The aim of this paper, as was stated, is to obtain the values of the different 
coefficients of the expression of energy due to the interaction between a solid wall 
and a liquid bulk. These values are associated with intermolecular potentials 
of the liquid and the solid wall. They take into account the molecular sizes 
and finite range of interactions between the molecules of the liquid and the 
solid. The value of coefficients 71 and 72 are positive. Consequently the 71 term 
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describes an attraction of the liquid by the solid, and the 72 term a reduction 
of the attractive interactions near the surface. In fact, the calculation leads to 
a more complex functional dependence than the one given by Cahn, Nakanishi 
and Fisher or do Gcnncs The energy of the wall takes into account not only 
the value of the density at the wall but also the gradient of density normal to 
the wall. This should be more accurate when the variation of density is strong 
with respect to molcciilar sizes. The density of the liquid i/i is chosen depending 
on three coordinates x, y, z. The direction to the wall is associated with z, and 
near the wall, the density may be chosen to be a function only of z. In fact 
such an assumption does not simplify the calculations. In our expression, the 
energy smoothly depends on x, y and allows continuous variations of density 
along the wall. The connection with Nakanishi-Fisher expression is not affected 
by this more general hypothesis. We denote by Es the new expression of the 
wall energy. 

The calculation may be extended to more complex cases: for example to the case 
of a liquid mixture in contact with a solid wall. We propose a general expression 
of the energy density at the solid surface. We note that the distribution of the 
concentrations of the components may be influenced by a solid wall effect. 



2. Description of the solid-fluid interaction in mean field approxi- 
mation 

In regions of a fluid where the density p in nonuniform, the intermolecular 
potentials produce a force on a given molecule that may generate surface tension 
effects Using classical molecular theory, it is possible to obtain a system 
of pressure and capillary tension that is mathematically equivalent to the stress 
tensor of a continuum model Such a description does not take into ac- 
count the possibility of interaction of the liquid with the solid wall (when the 
distribution of density varies strongly in the direction normal to a solid surface). 

These effects constitute the subject of the present paper. 

The effects are illustrated in Figure 1 and are described below. We consider a flat 
solid wall. The approach is aimed to describe the interaction of a macroscopic 
surface with a bulk fluid. The so-called "hard sphere diameter" of the molecules 
is denoted by a for the fluid and r for the solid. Then, the minimal distance 

between solid and liquid molecules is S = — (a + t) . 

In mean-field theory, we represent by V'(^) the intermolecular potential between 
two molecules of the fluid and respectively the potential between a fluid 

molecule and a wall molecule at separation r. The density of molecules at 
a point depends on its coordinates x, y, z, but the masses mi and TO2 
of each molecule of the liquid and solid are assumed to be fixed. The energy 
corresponding to the action of all molecules of the liquid and the solid on a given 
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molecule 1 located in (see Figure 1) is 

Wo = ^ ml il>iri)+Y^ 17111712 x{rj)- (1) 

i j 




Figure 1: Molecular layer between the liquid and the solid surface. 



The first summation (over i) is over fluid molecules (except for molecule 1), 

and the second summation (over j) is over the wall molecules. Molecule 1 is in 
the fluid and ri^rj are distances from molecule i or j to molecule 1. 

Denoting by viduji the number of molecules of fluid in the volume element 
dwi and i'2dw2 the number of molecules of solid in volume element dL02, the 
potential energy resulting from the action of all molecules in the medium on 
molecule 1 located in may be expressed in a continuous way: 

Wo = /// ml ipir)!/! dwi + /// mim2 x{r)j^2 dLJ2 (2) 
JJJni JJJn2 

where fli and CI2 are the domains occupied by the liquid and the solid. 

Potential energy Wq will be summed over all the molecules of the liquid. 
So, the first integral is counted twice in the previous integral over the domain 
Oi occupied by the liquid. In such a condition, we have to consider only the 
expression 

^ /// "^1 V'(?')f^i du-i + III mim2 x(r)j^2 di02 (3) 
^ JJJqi JJJn2 

for the potential energy Wq. 
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Let 1^1 (x, y, z) be an analytic function in each point M(x, y, z) of the Uquid. 

In the next derivation we replace at the point (0, 0, 0) Vi with its Taylor expan- 
sion in variables x, y, z limited to the second order. Such an assumption means 
that the forces between the solid and the liquid are of a short range. This is 
similar to the expansion given by Rocard^^. This means that il){r) is a rapidly 
decreasing function of r. Then, it is only necessary to know the distribution 
of molecules at a short distance from molecule 1. This case reflects the reality 
when the main force potentials decrease as with the distance. So, the po- 
tentials decrease very rapidly from the solid wall. The distribution of density 
inside the solid is assumed to be uniform. We write this expansion: 

v,{x,y,z) = -lo + x^ +,__+ (4) 

2 \ dx"^ ' dy"^ dz"^ dxdy dxdz dydz J 

where v\q , , ,■■ ■ represent the values of ui and its derivatives 

ux ox 
at point (0, 0, 0) . 

If we note pi = miVio and p2 = m2f2 the densities in the liquid at point 
(0, 0, 0) and inside the solid wall, we obtain 

/•oo poo 

Wq = 2tt rriipi / r'^tp{r)dr + n rriipi / r{h — r)tp{r)dr + 



r{h — r )tp{r)dr 



2tt rriipi / r'^tp{r)dr + n rriipi I 

J<7 Jh 

f°° r , ^ . X , dpi r 

2n mip2 / r{r - h)x{r)dr + - mi — 
Jh ^ c'z 

- miApi / r'^ip{r)dr H rniArpi / (-2r^ + 3hr^ - rh^)ip{r)dr + 

J mi ^ r r{h^ - r^)'4}{r)dr (5) 



In Eq. (5), h represents the distance between the molecule 1 and the solid 
wall, A is the Laplace operator and At is the Laplace operator tangential to 
the wall. 

Details of the calculation are given in Appendix 1. 

Notice that the last two terms may also be written in the form 

— miApi {-2r^ + 3hr^-rh^)ip{r)dr+-mi—^ / rh{h? -r^)il){r)dr 
12 Jh 4 dz Jh 

(6) 

The energy density per unit volume at point Q ]s Eo = vw Wq. 

In fact, this result is independant of the reference point. If we denote E the 
energy per unit volume at any point M in the liquid, we obtain 



5 



/■•oo />oo 

E = 27T pI r^iij{r)dr -\-7t pI / r{h - r)ip{r)dr -\- 

J <T J h 

roo ^ Opi 

27r P1P2 I r{r - h)x{r)dr -\- — pi / r{h? - r'^)ip{r)dr + 



h ^ Jh 

00 ^ roo 



TT _ 

Expression (7) yields 



/>oo 

i'^^^A (7) 



£ = i^p? + a(/i)pip2 + P{h)pl + jih)pi ^ + i^6ViApi + B (8) 
with 



/oo /"OO 

r^xl}{r)dr , Kb'^ = - J r'^^{r)dr 

POO POO 

a{h) = 27r / r(r — /i)x(r)dr , P{h) = tt / r(/i — r)'tp{r)dr 

Jh Jh 
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B=^piATPi {-2r^+3hr^-rh^)i;{r)dr+-p,^ r{h^-r^ 



)il){r)dr 



The constant 6 appearing in the piApi term denotes the covolume of the liquid 
as in the van der Waals equation ^^'^^ . 

The corresponding energy of the fluid is 

W = jjj EduJi (9) 



where E is given by Eq. (8) . 

To take into account the kinetic effects like in Rocard , p. 392, we must add 
the terms corresponding to kinetic pressure to the value of W. The first term 
Kp\ yields the internal pressure. 

Now, we calculate the different values of coefficients in Eq. 8 in a special case 
of London forces. 

3. Calculations of the energy of interaction in the case of London 
forces 

It is now possible to calculate the value of W for the particular interaction 
potentials. For example, one can take (see Appendix 2) 

^(0 = -i and x{r) = --^ (10) 
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In the case of London forces one has n=7. 

In fact this rough approximation is vahd only at short range from the wall 
(hypothesis ii). Following the calculations in Appendix 2, we obtain that the 

two integrals in the B term are of the order of — ; moreover, 



a{h) 



V2h? 



h 



(11) 



Integrals associated with Eq. (9) are taken over an interval [S, d\ , where d is 
the range of molecular forces in the liquid and S is the minimal distance between 
solid and liquid molecules. 

For the same reasons as in Eq. (4), the expansion of pi , , at the 

oh ah-^ 

solid wall is taken in the form 



Pi = 

dpi 
dh 

dh? 



Ps 



dp 



dps 
dh 

d-'ps 
dK^ 



: h^ 
dh ^ 2 
d'ps 



d'ps 

dh^ 



+ 



+ h 



+ 



dh'^ 



+ 



(12) 



Terms Ps , ^ denote values of pi and its normal derivatives at the 

dh dh'' 

solid wall. 

rp,. / N dps ,dpi d'^ps fd^Pi^ 

This means ps = {pi)s, ^ = i^)s, ^ = i^)s. 

Now we make the important assumption that the variations of pi take into 

5 

account several molecular ranges. Hence, - can be considered as a small pa- 
rameter. It implies that the first derivative of pi with respect to h is on the 

order of — , the second derivative of pi is on the order of Then, the B 
a 

term can be removed from the integration of E. Consequently, in the calcu- 
lation of the surface energy in Eq. (4), the second derivative terms may be 
removed (but not in the calculation of the bulk energy of the liquid associated 
with K{pI + b^piApi)). This result is in agreement with some considerations 
by de Gennes 

Keeping terms up to the first order in Eq. (8) yields 



E = K{pl+b'piApi)+a{h)p2 +P{h) (^pl+2hps^^ +^^^^P'^ 



or 



E = K{pl+b'^piApi)+a{h)p2Ps+P{h)ps+ha{h)p2 



dps ^ 2hp{h)+j{h) dpi 



dh 



+ 



dh 
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Then, the energy of the fluid is 

W = Wi+W2 , with 
Wi = fff K{pl + b''piAp,)du;i 



W2 = P2PS III a{h)dhdxdy + Ps /// l3{h)dhdxdy 



with 



a{h)dh = , J^^°° mdh ■■ 



k-K 
2AP 



k-K 
V25 



j{h) kn 



Wc note that piApi = pi div{grad pi) = div{pi grad pi) — [grad pi) 
and the Stokes formula yields 

^ IIL ^^''^ ~ ^'^^ Pi?)d^i + jj^ Kb'' PS ^ ds (14) 

Here S notes the surface of the solid wall corresponding to the boundary of 
fii. As it was said in paragraph 2, to take into account the kinetic effects, we 
must add to energy Wi the terms corresponding to kinetic pressure. 

The first term Kp\ yields the internal pressure, and consequently 

e = a{pi,r]) - {grak pif 

Pi 

denotes the internal specific energy. Term pia{pi,ri) is the bulk internal volume 
energy as a function of pi and of the specific entropy rj in the liquid. 

Let us note that 

Then, a straighforward calculation yields the energy of the liquid in the well 
known form 
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but with 



dp% „ dps 

Here, 74 means 2^574 and 



71 = P2, 72 = ^ (16) 

/LtTT kir f 1 1 



73 = ^ P2, 74 = g ^ ^ g ^ 

is the form of the special energy to be added at the sohd surface to obtain 
the total energy of the liquid. In expression (15), the 71 term (favoring large 
ps) describes an attraction of the liquid by the solid. The 72 term represents a 
reduction of the liquid/liquid attractive interactions near the surface : a liquid 
molecule lying directly on the solid does not have the same number of neighbors 
that it would have in the bulk. The terms with the coefficients 73 and 74 also 
describe a reduction of the liquid/liquid attractive interactions due to the lack 

d 

of molecules of the liquid near the wall (in the case where — — is positive). 

oh 

Expressions (15) and (16) do generalize the results by Nakanishi and Fisher^: 

Expression (15) contains terms (the first two) similar to the expression given by 
Nakanishi and Fisher, but we find additional terms associated with 73 and 74. 
Additional terms are the correction of the two first terms. Our main interest 
is to obtain values of coefficients of the energy of the wall as a function of the 
properties of molecules and to take into account variations of density at the wall 
strong enough with respect to molecular sizes. 

4. Extension to the case of liquid mixtures in contact with a solid 
wall 

Here we propose an extension of the above theory to the case of liquid mixtures. 
This example is for a binary mixture, but there is no reason one could not include 
more species in the mixture. 

The hypotheses are the same as in the case of one-component liquid. The 
only difference will come from the interaction of molecules of the two liquids. 

An adaptation of the previous calculations yields the following results. The 
potential energy resulting from the action of all molecules in the medium on 
molecule 1 of liquid 1 located in is 

WiQ = /// ml ■>pi{r)vi(kji 
JJJni 



+ 



/// mim2ips{r)v2dhJi + /// mimsXi{r)i'sdhJ2 (17) 
JJjQi JJJq.2 
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This energy is only for one species in the Hquid. To determine the whole energy, 
one must first sum over all molecules of species 1 and then do an analogous 
summation over the molecules of species 2. Potential energy H^io will be 
summed over all molecules of the liquid mixture. In this way the contribution 
of the liquid-liquid integral, Wio, will be counted twice in the previous integral 
over the domain Cli. 

We have denoted by m^, the mass of the molecule of fluid i (i e {1,2} ), tpi 
and V2 are the potentials of interaction between molecules of fluid 1 with 
themselves and among molecules of fluid 2, is the potential of molecular 
interaction between the two fluids, Xi (* G {li^}) are the intermolecular po- 
tentials of fluids with the solid wall, and (i € {1, 2, 3} ) denote the number of 
molecules in liquids and solid wall per unit volume. 

Then, following the procedure of section 2, we obtain : 

/•OO 

Elo = z^ioH^io = 27r Pi / r'^'4}i{r)dr 

J (Tl 

/•OO fOC 

+ -K p\ I r{h — r) ijji{r)dr + 2iTpipz / — h)xi{r)dr 

Jh Jh 

+ ^ Pi ^ J r{h? - r^) il}i{r)dr + ^ pi Api j r^tpi{r)dr 
+ — Pi AtPi / (-2r^ + 3hr^ - rh^)il)i{r)dr 

/■OO POO 

+ 2 'Kpip2 I r'^ips{r)dr + 7rpip2 / r{h - r)i)^{r)dr 

TT d poo ^ poo 

+ I Pi r{h'^ - r^) ip3{r)dr + ^ Pi ^P2 J r'^'^i{r)dr 



f-OC 

/ (- 
Jh 



+^ pi AtP2 / {-2r^ + 3hr^ - rh^)^i{r)dr 



where cr,; (« £ {1,2}) is the diameter of molecule of liquid i. 

We may also repeat the same calculation for £^20 associated with the second 
component of the mixture. As for a simple fluid, we may take into account 
the kinetic effects and the internal energy of a nonhomogeneous mixture as in 
Fleming et al Then, we obtain the following additional energy at the solid 
surface for the liquid mixture in the form 

Es = -711 Pis - 721 P2S + ^ ( 712 Pis + 722 P2S + 2732 PISP2S ) 
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_^,3 ^ -723 ^ - f 714 ^ + 724^ + 2 734 ^^^^ (18) 

dh dh \ dh '^^ dh dh J ^ ' 

All the coefficients 7^^ can be calculated explicitly after the particular form 
of interaction potentials was chosen. For example in the case of London forces, 
the values of coefficients related to the densities of the two fluids at the surface 



are 



= 125? = 12lf 

= 125?' = 12lf 



where ps is the density of the solid, /Xj are the coefficients associated with po- 
tentials Xi) i € {1,2}, ki is associated with potentials tpi, i e {1,2,3} and 

Si = - (cti + r), i e {1, 2} are the minimal distances between the solid and 

molecules of the two species of the mixture. 

Such an expression allows one to estimate the influence of a solid wall on each 
component of a fluid mixture. Depending on the values and signs of different 
coefficients -jij , one can estimate the magnitude of the attraction or repulsion 
effects due to the wall. Placing the fluid mixture in contact with specially 
chosen solid walls may be an efficient way to separate constituents of a molecular 
mixture. 



5. Conclusion 

The energy of a fluid in contact with a solid wall contains a contribution 
from the solid which may be represented by a surface density function. For a 
flat wall wc obtain this expression by taking into account not only the density 
of the liquid at the surface but also its normal derivative. All the numerical 
values of the coefficients of the surface energy functional are calculated in terms 
of the parameters of molecules in solid and liquid. This energy characterizes the 
behavior of the surface in contact with the fluid in the wetting transition. The 
method may be extended to the case of nonflat solid surfaces which is important 
in catalysis chemistry. 

Appendix 1. 
Calculation of the value of E in Eq. (8) 



To obtain the formula for E given in Eq. (8), we have to calculate the two 
integrals: 
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^ III aiid JJJ mim2x{r)i^2dLi)2 

a) Calculation of ^ J^j^ ■;/'('') ^^ic^^^i 

This integral is associated with the energy of interaction between molecules 
in the liquid. Let us denote by S{0,h), the domain occupied by the sphere 
centered at (0, 0, 0) and with radius h (see Figure 1). We introduce {r,6,(fi) 
the spherical coordinates associated with the center of the sphere. Then, 

^ Ilia ™^ "^^^^ " 
t; Iff iTT-l V(^) ciwi + i Iff ml il){r) ui dui 

Let us note than for any integers p, g, r, and any boundary S of a sphere 
centered at (0, 0, 0) and radius r, 



47r 4 



// x^f+i y« ^'^ ds = and // ds = [[ ds = [[ ds = 
JJi: JJn JJn JJt. 

Then, 

\ 1(1 m\'i()(r)vidLOi = 2'iTrn\viQ f r'^tlj{r)dr + ^mlAvio ( r'^ip(r)dr 
2 JjJs{o,h) J a 3 

and 

rn{il}{r)vidw\ = - I i j ( j m1ip{r)i'isiinpd(fi\r^dr 

2 Jo \ Jh ^ JArccos!^ ' I 



r-2-n I r- + oo 

l(li~S{0,h) " ' Jo \Jh ^JAtccos^ 

The expansion up to the second order of Vi (see Eq. (4)), yields 



^ /// ™i ^('') '^'^i = IT uio ml [ r {h + 

^ JJJai-S{0,h) Jh 

POO 

I r{h^- r^) iPir) dr 

Jh 



r) 'tjj{r) dr 



TT dl^io 2 

+ 2 



TT / d^Ulo , d^l^io 1 2 

— ^ + I " 

12 \ dx^ dy 



+ ^ I + ] „j2 (2r^ + 3r^ h-rh^) i,{r) dr 



+ l^rnl r{h' + r')^{r)dr 



which is the desired result. 
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Figure 2: Representation of the changed variable associated with the molecular 
layer at the solid surface. 



b) Calculation of / / / mi m2 x(^) ^2 duj2 

JJJQ.2 

This term corresponds to the energy with respect to the solid wall. 



p27V / /-CXD 



mi m2 xir) ^2 dw2 

'Arccos ^ 



= mi m2 1^2 J yj x('') smip d(pjr'^ dr \ d6 

POO 

= 2-K mi m2 f 2 / r {r — h) x(r) dr 

Jh 

which gives relation (7). 

Now, we change variable as in Figure 2. The origin of the third coordinate z 
is placed at 0'. Consequently z = O'O with the position of molecule 1. 

If wc take now the origin of the z-axis at the solid wall and change the 
orientation in such a way that 00" = h, we obtain that z + h is a constant 
for a given molecule. Hence, we get 



Then, 



dpi _ _ dpi d'^pi _ d'^pi 
dh dz ' dh'^ dz^ 



POO POO 

E = 2-Kpl / ip{r) dr + npl / r{h - r) 4>{r) 

J a Jh 
TT dpi 



dr 



1*00 ^ Q I'OO 

+27rpip2 J r{r - h) x{r) dr + - pi J r(h? - r'^) i){r) dr 
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(we note similarly the Laplace operator in new coordinates and in old coordi- 
nates). So, we obtain Eq. (8). 



Appendix 2. 

Some remarks about the potential associated with cohesive forces 

In molecular theory it is proved by the virial method than the so-called "coeffi- 
cient a" of the van der Waals equation can be obtained from 



f 

J a 



fir) Anr^ dr = (20) 



The term f{r) is the magnitude of the attractive force between the two 
molecules of the fluid at the distance r, and TV is the Avogadro number. The 

potential V'i(^) ^ ^ J fi^) dr is assumed to be zero at infinity. 

Take the forces in the form 

= 4 ' Mr) = ^ 



r" ' ^ ^ (n - l)r"-i 
then, 



A A r 3-n ^ 47r A 1 
Aw A r"^ " dr 



I 

J a 



n - 4 iT"-4 
andEq. (20) yields 

_ 3(n-4)aa"-^ 

^ - 2 TT iV2 i^l) 

Let us introduce tpir) such that ml ip{r) = ipi{r). The term ip gives the 
potential per unit mass. Then, 

-3(n-4) a cr"-^ 

tp{r) = 



2(n -1) TT ml r"-i 

Here, Nm = M denotes the molar mass of the fluid. In case n = 7 (London 
forces), we get: 

' A-KM^r^ 

In the following, we use: 

, , , k . , , 3a (7^ 

rv) = — fi with k 



An M2 
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Consequence: Calculation of E in Eg. (8) for London forces. 



Now, using direct integration, we obtain the values of the coefficients given 
in Eq. (8). Since now 



we obtain 



POO 

27r / tp{r) dr 

J a 



«n— 1 



2 k 7 
"3^ 



and n = 1 



ip(r) dr = — - — 



k IT 
3~a 



h) x(r) dr 



2 TT / r (r 

Jh 

poo 

■K r {h — r) tp{r) dr 
Jh 



jJL TT 



k TT 

127^ 



f°° /, 2 2x , / X , ^ 



TT 1'°° k TT 

— / (-2r^ + 3hr^-r h^) Mr) dr = —- 
12 ift ^ ' 16 h 



r {h'^ - r"^) ■>p{r) dr 



k TT 
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